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We examine the lowest excitations of one-dimensional few-boson systems trapped in double wells of variable 
barrier height. Based on a numerically exact multi-configurational method, we follow the whole pathway from 
the non-interacting to the fermionization limit. It is shown how, in a purely harmonic trap, the initially equidis- 
tant, degenerate levels are split up due to interactions, but merge again for strong enough coupling. In a double 
well, the low-lying spectrum is largely rearranged in the course of fermionization, exhibiting level adhesion and 
(anti-)crossings. The evolution of the underlying states is explained in analogy to the ground-state behavior. 
Our discussion is complemented by illuminating the crossover from a single to a double well. 
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I. INTRODUCTION 

The realization of Bose-Einstein condensates has made ul- 
tracold atoms an ideal tool forprobing and understanding 
paradigm quantum phenomena IHIJjHH]- One such example 
is the quasi-one-dimensional Bose gas, where the transverse 
degrees of freedom are frozen out such that an effective one- 
dimensional description becomes possible. As it turns out, in 
such a system one can tune the effective atom-atom interac- 
tion strength at will by merely changing the transverse con- 
finement length |5]. This allows us to explore the limit of 
strong correlations and, in this way, to test the physics beyond 
the mean-field description applicable to large and weakly in- 
teracting systems. 

The limit of infinite repulsion — so-called hard-core 
bosons — is in turn known exactly in one dimension, as it is 
isomorphic to that of free fermions (up to permutation sym- 
metry la]). This makes it appealing to think of the exclusion 
principle as mimicking, as it were, the hard-core repulsion, 
which is why this limit is termed fermionization. Its theo- 
retical description has thus attracted a great deal of attention 
1 3, 1 MB . [Toll , stimulated by its recent experimental verification 

inHl. 

On the other hand, the question how exactly these two very 
different borderline cases connect has attracted less attention. 
The major reason for this is that it is notoriously hard to in- 
clude the effects of strong coiTelations from first principles. 
This can only be done forrather few particles, say N ^ 10. As 
it happens, this is not only the regime of experimental acces- 
sibility 1 13]; it is also very interesting because the signatures 
of two-body interactions are still pronounced, thus facilitating 
the understanding of larger systems. An analytic solution is 
known for homogeneous systems with periodic |3 HH] and 
fixed boundaries [16], but also for the elementary case of two 
atoms in a harmonic trap |17], which already captures some 



key features of the evolution from the weakly to the strongly 
interacting limit. In general, though, the solution of trapped 
interacting bosons requires numerical approaches. Most in- 
vestigations have so far focused primarily on the ground state 
ifisi Il9[ I20I I211 112, HU, or on regimes where correlations are 
weak enough to be p assably represented by few single-particle 
orbitals Im [111261. 

The goal of this paper now is to extend the systematic study 
of the fermionization transition in harmonic and double-well 
traps to the lowest excitations of finite boson systems. Un- 
derstanding the low-lying spectral properties is not only an 
interesting problem in its own right, given the richness of the 
pathway to fermionization for the ground state. It is also an 
essential contribution to the control of few-body systems (as 
is desirable for quantum information processing) and to ex- 
plaining their dynamics, where the double well is a prototype 
model for fundamental effects such as tunnehng or interfer- 
ence '28', '2^ . As in our previous works 23], we draw 
on the Multi-Configuration Time-Dependent Hartree method 
lHOillllll^]- As its single-particle basis is variationally opti- 
mized for each state, it allows us to study the excited states of 
few bosons in a numerically exact way. 

Our paper is organized as follows. Section HIl introduces 
the model and gives an overview of some key concepts. In 
Sec. Hn] we give a brief introduction to the computational 
method and how it can be applied to the study of excitations. 
The subsequent section finally is devoted to our results for 
the low-lying spectrum (Sec. lIV Al l and the underlying excited 
states (Sec. lIV ST i. Our discussion is rounded off by illuminat- 
ing the crossover from a single well to a double-well trap in 
Sec. UVCl 

II. THEORETICAL BACKGROUND 
A. Model 
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In this work we investigate a system of few interacting 
bosons {N — 2, . . . , 5) in an external trap. These particles, 
representing atoms, are taken to be one-dimensional (ID). 
More precisely, after integrating out the transverse degrees of 
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freedom and upon introducing dimensionless variables we ar- 
rive at the model Hamiltonian (see [22] for details) 



where h = + U{x 



is the one-particle Hamiltonian with 
a trapping potential U, while V is the effective two-particle 
interaction potential iH] 



V{x) = gSaix), with 5 
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Here an s-wave scattering length gq and a harmonic transverse 
trap with oscillator length a± were assumed. The well-known 
numerical difficulties due to the spurious short-range behavior 
of the standard delta-function potential S{x) are alleviated by 
mollifying it with the normalized Gaussian 



S^{x) 
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which tends to 6{x) as cr in the distribution sense. We 
choose a fixed value a = 0.05 as a trade-off between smooth- 
ness and a range that is much shorter than the length scale of 
the trap, oii — 1. 



B. Correlations and fermionization: key aspects 

Our approach equips us with the full solution of the 
system — ^here, the excited-state wave functions, which are 
fairly complex entities. Visualizing and in this way relat- 
ing them to the physical picture, it is useful to consider re- 
duced densities, or correlation functions. As is well-known, 
the knowledge of some wave function is equivalent to that 
of the density matrix = To the extent that we 

study at most two-body correlations, it already suffices to con- 
sider the reduced two-particle density operator 



P2 ^ tr3..Ar|*)(*| 



(1) 



whose diagonal kernel P2{xi,X2) gives the probability den- 
sity for finding one particle located at xi and any second one 
at X2- For any one-particle operator, of course, it would be 
enough to know the one-particle density matrix pi = tr2P2, 
so that the exact energy may be written as 

E = Ntr{p,h) + ^^i^l—2ltr{p2V). 

The one-particle density matrix can be characterized by its 
spectral decomposition 



Pi 



= ^na\(l)a){<l)a\ 



(2) 



Remarkably enough, not only the non-interacting limit is 
well known, but also the complementary case of infinitely 
strong correlations, 5 ^ oo. It is commonly referred to as the 
Tonks-Girardeau limit of ID hard-core bosons, or also as their 
fermionization. This lingo finds its justification in the Bose- 
Fermi map [33ll that establishes an isomorphy between the 
exact bosonic wave function and that of a (spin-polarized) 
non-interacting/ermion/c solution 'I'g^, 
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(3) 



where A = Wi^j sgx\{xi — Xj). The mapping holds not only 
for the ground state, but also for excited and time-dependent 
states. Since = 1, their (diagonal) densities as well as 
their energy E will coincide with those of the corresponding 
free fermionic states. That makes it tempting to think of the 
exclusion principle as mimicking the interaction (g 00), as 
is nicely illustrated on the ground state of N hard-core bosons 
in a harmonic trap |9] 



vl/+(g)(xe- 



where Q ^ {xi, . . . ,xn)^ ■ 



l<i<j<N 



III. COMPUTATIONAL METHOD 

Our goal is to investigate the lowest excited states of 
the system introduced in Sec. |II] for all relevant interac- 
tion strengths in a numerically exact, i.e., controllable fash- 
ion. This is a highly challenging and time-consuming task, 
and only few such studies on ultracold atoms exist even for 
model systems (see, e.g., |l2l1, ^2^). Our approach relies on 
the Multi-Configuration Time-Dependent Hartree (mctdh) 
method |32], primarily a wave-packet dynamics tool known 
for its outstanding efficiency in high-dimensional applica- 
tions. To be self-contained, we will provide a concise in- 
troduction to this method and how it can be adapted to our 
purposes. 



A. Principal idea 

The underlying idea of MCTDH is to solve the time- 
dependent Schrodinger equation 



as an initial-value problem by expansion in terms of direct (or 
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Hartree) products = ipj_^' 

vI/(Q,i) = ^A7(t)<i>j(Q,t). 



(4) 



where Ua £ [0, 1] is said to be the population of the natural or- 
bital (pa- If the system is in a number state |n) = \nQ,n[, . . .) 
based on the one-particle basis {(pa}, then = n'a/N in (|2]); 
but it also extends that concept to non-integer values of nj,. 



The (unknown) single-particle functions ip^^^ (j^ — 
1 , . . . , n^) are in turn represented in a fixed primitive basis 
implemented on a grid. For indistinguishable particles as in 
our case, the sets of single-particle functions for each degree 
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of freedom k = 1, . . . , iV are of course identical (i.e., we have 
(fj^, with j„ < n). 

Note that in the above expansion, not only the coefficients 
Aj are time-dependent, but so are the Hartree products $j. 
Using the Dirac-Frenkel variational principle, one can derive 
equations of motion for both Aj,ipj ll32ll . Integrating this 
differential-equation system allows one to obtain the time evo- 
lution of the system via (|4]i. Let us emphasize that the concep- 
tual complication above offers an enormous advantage: the 
basis {$ J (•,<)} is variationally optimal at each time t. Thus 
it can be kept fairly small, rendering the procedure very effi- 
cient. 

It goes without saying that the basis vectors $ j are not per- 
mutation symmetric, as would be an obvious demand when 
dealing with bosons. This is not a conceptual problem, 
though, because the symmetry may as well be enforced on 
^ by symmetrizing the coefficients Aj. 



B. Application of the metliod 

The MCTDH approach 134] . which we use, incorporates a 
significant extension to the basic concept outlined so far The 
so-called relaxation method fss'] provides a way to not only 
propagate a wave packet, but also to obtain the lowest eigen- 
states of the system, '^m- The key idea is to propagate some 
wave function ^'^^^ by the non-unitary e~^'^ (propagation in 
imaginary time.) As r ^ oo, this exponentially damps out 
any contribution but that stemming from the true ground state 
like e~^"^. In practice, one relies on a more sophisticated 
scheme termed improved relaxation |30, 36], which is much 
more viable especially for excitations. Here {'^\H\'^) is mini- 
mized with respect to both the coefficients Aj and the orbitals 
ipj. This leads to (i) a self-consistent eigenvalue problem for 
which yields Aj as 'eigenvectors' , and (ii) 
equations of motion for the orbitals Lpj, but based on certain 
mean-field Hamiltonians. These are solved iteratively by first 
diagonalizing for Aj fixed orbitals and then 'optimizing' 
ifj by propagating them in imaginary time over a short period. 
That cycle will then be repeated. 

Whereas the convergence to the ground state is practi- 
cally bulletproof, matters are known to get trickier for excited 
states (see 13611 ). This should come as no surprise, granted 
that one cannot just seek the energetically lowest state possi- 
ble but should remain orthogonal to any neighboring vectors 
(That is why, at bottom, the convergence turns out to 
be highly sensitive to the basis size — that is, to n — even for 
small correlations: The lower states simply must be repre- 
sented accurately enough.) For practical purposes, the most 
solid procedure has proven to be the following. In the non- 
interacting case, we construct the eigenstates as number states 
\n) = |nQ, n']^, . . . ) in the single-particle basis {(t>a]- Starting 
from a given |n), the eigenstate 5',„ for g > is found by an 
improved relaxation while sieving out the eigenvector closest 
to its initial state \n). The resulting eigenstate will then in turn 
serve as a starting point for an even larger g, and so on. 

As it stands, the effort of this method scales exponentially 
with the number of degrees of freedom, . This restricts our 



analysis in the current setup to about = 0(10), depending 
on how decisive correlation effects are. Since the computation 
of excited states requires that the neighboring states be suffi- 
ciently well represented, the basis must in fact be rather large 
even for weak correlations. 

As an illustration, we consider systems with 5 and 

need n ~ 15 orbitals. By contrast, the dependence on the 
primitive basis, and thus on the grid points, is not as severe. 
In our case, the grid spacing should of course be small enough 
to sample the interaction potential. We consider a discrete 
variable representation with 95 to 125 grid points per degree 
of freedom. 



IV. LOWEST EXCITATIONS 

As in Refs. ||22[ l23ll . we consider bosons in a double-well 
trap modeled by 

U{x) = ia;^ + h5w{x), 

expressed in terms of the harmonic-oscillator length a\\. This 
potential is a superposition of a harmonic oscillator (HO), 
which it equals asymptotically, and a central barrier which 
splits the trap into two fragments. The barrier is shaped as 
a normalized Gaussian S^j of width w and barrier strength 
h. If w ^ 0, the effect of the barrier reduces to that of a 
mere boundary condition (since 5w S), and the correspond- 
ing one-particle problem can be solved analytically |37]. Al- 
though this soluble borderline case presents a neat toy model, 
the exact width w does not play a decisive role and is set to 
w = 0.5. 

In Refs. ll22l l23ll . we have studied the ground-state evo- 
lution from the weakly correlated regime to fermionization, 
with an eye toward the fascinating interplay between inter- 
atomic and external forces as the barrier height h was varied. 
We now seek to extend that investigation to the lowest excita- 
tions. In Sec. IIV Al we will look into the low-lying spectrum 
a{H) = {E„i}, whose corresponding eigenstates '^„i will be 
analyzed in detail (Sec. lIVBb . As the spectral properties in 
the cases of a single and a double well will turn out to be quite 
different, the question as to how they connect naturally arises. 
That crossover will be the subject of llV CI 

A. Spectrum 

In this section, we study the evolution of the lowest energies 
Em (g) as g passes from the non-interacting to the fermioniza- 
tion limit. Figures [Tl3] convey an impression of this transition 
for N = 3, 4, 5 bosons in a harmonic trap (h = 0) and in a 
double well (h = 5), respectively. Before dwelling on the de- 
tails, let us first capture some universal features of the spectra. 

In the uncorrected limit, g — > 0, the energies are simply 
given by distributing the atoms over the single-particle levels 
Eq, starting from tiq — N (the Bose 'condensate'): 

S = iVtr(pi/i) = ^<e,. (5) 

a 
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In particular, _Bo = iVeo; hence the 'chemical potential' I-In = 
= eo, as usual. Note that Eq. ^ implies 



^0 - ^0 



degeneracy if two single-particle energies are commensurate, 
i.e., J^ai'^'a - ^a)eQ = for two n^fi. 

In the Tonks-Girardeau limit, on the other hand, the spec- 
trum becomes that of a free fermionic system (even though, of 
course, the system is really still bosonic and has an all but neg- 
ligible share of interaction energy). Thus one can find some 
(auxiliary) n with e {0,1} such that 



lim = > 7i; 



(6) 



In the ground state, the particles can therefore be thought of as 
filling the energy ladder up to the Fermi edge, ea < £n = fJ-N- 
For a harmonic confinement, the chemical potential will thus 
be cx N, so = 0{N^). 

It should be pointed out that, in the spirit of the Bose-Fermi 
map (O, the borderline cases of no and infinite repulsion may 
be perceived as one and the same (non-interacting) system, 
their sole difference being the 'exchange symmetry' emulat- 
ing the effect of interactions. Therefore the same type of en- 
ergy spacings and (quasi-)degeneracies should appear at both 
ends of the spectrum. 
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Harmonic trap (h — 0) 

For a single well, the one-particle spectrum {ca — a + ^} is 
known analytically, which readily equips us with the full spec- 
trum for both the non-interacting and the fermionization limit. 
First consider the case g = 0. Then Eq — N/2, while all 
other levels follow with an equal spacing of Aq = 1. Owing 
to that equidistance, the degree of degeneracy # = Em — Eq 
(m > 1) goes up with each step, measured by the average 
occupation Na = J2a "-'a'^- Explicitly, while both m — 0,1 
are non-degenerate, the eigenspace pertaining to E2 — E3 — 
N/2 + 2 is two-dimensional (see Fig.[T]i, etc. 

To understand this degeneracy and how it is lifted, let us 
recall that, in a harmonic trap with homogeneous interactions 
V{xi — Xj), the center of mass (CM) R :— X^iLi ^i/^ is 
separable from the relative motion. Hence one can decompose 
the Hilbert space H = Hcm ® Hid so as to write 

^ V^rci; E{g) = (AA + i) + e,^i{g). 

This signifies that for every level for the relative motion, 
e^-ciig), there is a countable set of copies shifted upward by 
A/" = 0, 1, . . . . For g — 0, V-'rci is a harmonic eigenstate as 
well, so e^^j^ {Q) + for some V, and several different 
combinations of {M , v) may coincide. Switching on g > 0, 
however, breaks that symmetry, leaving M untouched while 
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Figure 1: Lowest energies Em in a harmonic trap {h — 0) for A'^ = 
3, 4, 5 bosons. (The lines connect the data points to guide the eye.) 



This fact is nicely illustrated on the example of = 2 atoms, 
where S 



dg 



^rcl 



= (V'.|<5(r)|^,) = 1^,(0)1' CX 



( ) denoting the binomial coefficient; so higher excited rela- 
tive states 'feel' the interaction less. This fits in with our find- 
ings in Fig. [T] The two states m = 2, 3 break up, the lower 
curve — in light of the reasoning above — ^pertaining to higher 
internal excitation. 

Apart from that, the spectral pattern does not give an air 
of being overly intricate but follows the general theme known 
from the two-atom case. All levels first rise quickly in the 
linear perturbative regime, but start saturating once they enter 
the strongly interacting domain {g ^ 10). As insinuated, the 
fermionization limit is known exactly, which endows us with 
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Figure 2: (color online) Single-particle spectrum {ea} in a double 
well with barrier height h = 5. 



a helpful calibration. Since the limits g 0(oo) can be re- 
garded simply as bosonic (fermionic) counterparts of the same 
non-interacting system, the two should share exactly the same 
energy scales (here Aq — 1). Indeed, building on the ground- 
state energy Eq — J2a<N ^i ~ all levels again follow 
in equal steps Aq, with a degeneracy # = Em — Eq. This 
fact, effortless as it may come out of the theory, is a strong 
statement, for it implies that the very interaction that divorces 
some degenerate lines at g = is also responsible for glu- 
ing them together again if it gets sufficiently repulsive. An 
indication of this effect may actually be observed in Fig.[T] 



Double well (h = 5) 

As opposed to the purely harmonic trap, the single-particle 
spectrum {ea} of the double well is not that simple. In order 
to get a rough idea, it is legitimate to consider a toy model 
of a delta-type baiTier h6{x) (i.e., w 0) |37]. Then only 
the even single-particle functions (f)a will be affected as they 
have nonzero amplitudes at a: = 0. For h > 0, these will 
be notched at zero, and in the limit of large enough baiTiers, 
h ^ oo, their density will approach that of the next (odd) 
orbital (pa+i — while still remaining even — and so will their 
energy, Ca fa+i from below. In that extreme case, we 
would end up with a set of doublets separated by gaps Ae = 2. 
The non-interacting many-body spectrum {En = X^a'^a^o} 
would then be composed of a lowest band of states within the 
(1 + A^)-dimensional subspace span{|nQ, n'j^ = N — n'^) \ 
Uq = 0, . . . N}, while the next band (obtained by removing 
one particle from the lowest levels eo/i) would be shifted up- 
ward by Ae. 

The realistic single-particle spectrum is sketched in Fig.E] 
Due to the nonzero width w, also the odd states are shifted 
slightly, and of course the distance between the doublets is 
never exactly Ae = 2 as in the crude estimate above. What 
is more, the finite barrier h = 5 not only lifts the even-odd 
degeneracy already for the lowest states, where the doublet 
character is still pronounced and the gap A/, = ei — eg ^ 
1, but even more so for levels above the barrier. There the 
central region is classically allowed and the spectrum becomes 
more and more harmonic for higher ea- The consequences for 
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Figure 3: Lowest energies Em in a double well (h — 5) for TV = 
3, 4, 5 bosons. Inset (A'' = 4): level adhesion for the states m = 3, 4 
(counted from below al g = 0). 



the full non-interacting spectrum a {H) in Fig. [3] are that the 
lowest band is 'fanned up' into {E,n — Neq + mA/j}. The 
next cluster of levels is still well separated in energy. 

The situation gets slightly more involved in the fermion- 
ization limit g ^ oo. Here the spectrum is generated by 
(fictitious) fermionic states \n) with G {0, 1}, so En = 
^'a^a is obtained by removing particles from below the 
Fermi edge en — fJ-N- The first excited state would thus 
be higher in energy by eN — £n-i, which — in our simplis- 
tic model invoked before — would be about zero for odd N 
and two for even N. The next band would then be made up of 
four levels for any number of atoms, created by pushing one 
particle out of the doublet pertaining to e^r-i into the next 
upper doublet. However, since these involve higher energies. 
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the typical doublet structure encountered in the lowest levels 
is lost for our finite barrier h = 5, and the spectrum will be 
rearranged and mixed, as seen in Fig. [3] Here the bands are 
no longer separated but smeared out considerably. Still, the 
odd-vs-even distinction with respect to the lowest level may 
be observed. 

How do these two ends of the spectrum connect? As can be 
seen in Fig.[3j the reordering that is necessitated by the Bose- 
Fermi map starts taking place already for small g, where the 
influence of the interaction operator Hi := J2i<j ^{^i ~ ^j) 
can still be treated as a perturbation. 



AS„(.g) ^ {Hi) 



^ \{n'\Hi\n)f 
^ ^„(0)-S„,(0)- 



In light of this, the different effects of Hi on different (non- 
interacting) states \n) will manifest themselves in (i) their 
slopes and (ii) different curvatures. While the slopes £"(0) > 
do not differ dramatically, the second order fosters level re- 
pulsion within each band, since the energy gaps are small in 
magnitude (|£'„(0) (0)1 0{Ah=5) < 1) but have 
different signs. (Needless to say, by conservation of parity 
n^((5) = Q) solely states within the same symmetry 
subspace {^I' | 11^1' = ±5'} are coupled.) Even though lim- 
ited to perturbative values g < 1, this may help us under- 
stand why the upper lines of the lowest band tend to rise so 
steeply, whereas the lowest ones in each band are 'pushed' 
downward a little. For example, see Fig. ^ N ~ 4: the two 
levels m — 3,4 (counted at g — 0) intersect the lowest upper- 
band state |3, 0, 1) at g ~ 1. 

Most stunning is the observation that, apparently, some 
lines are virtually 'glued together' once the interaction is 
switched on (see Fig. [3j insets iV = 4, 5). In particular, this 
applies to the upper pair m ^ N — 1,N (counted at g ~ 0, 
as always). To be sure, the two levels are close from the 
start (AE = A5 ~ 0.1); but for g ^ 1 they get as close 
as ^ 0.01. This quasi-degeneracy arises as interactions are 
turned on; but it is destroyed again once these get very strong 
(g ~ 5). However, there is no indication as to what exact 
mechanism brings the two lines together Not only do the cor- 
responding states remain orthogonal at all couplings, but they 
stem from opposite-parity subspaces and can thus only mix 
with other states of a kind. Still, their reduced densities will 
turn out very much alike (as will be laid out in Sec. lIV Bl i. 



B. Excited states 

As yet, we have looked into the spectrum and its evolu- 
tion from the weakly to the strongly interacting regime. In 
Sec. IIV A| we have found characteristic spectral patterns for 
the cases of a single and a double well, respectively. We now 
aspire to get a deeper insight into the underlying states 5'm>i, 
which may be also beneficial for studying the dynamics in fu- 
ture applications. 

Generally speaking, the non-interacting limit is described in 
terms of number states \n) in the respective one-particle basis. 
Owing to the asymptotically harmonic confinement, we thus 



have an overall Gaussian profile p{x) oc exp (— a;^), which 
is modulated by the central barrier as well as the degree of 
excitation. At least for the low-lying states, the length scale 
is therefore about that of the harmonic confinement, a|| = 1. 
Being single-particle states, they are essentially devoid of two- 
body correlations, which amount to /92 = ^(1 + Pi2)pi ^ pi 
(with the permutation operator P12). 

When interactions are added, some extra interaction energy 
^^^'"'"'' tr {Vp2) must be paid. Hence, the system will re- 
spond by depleting the correlation diagonal P2{xi, X2 — xi), 
roughly speaking. In our single-particle expansion, this will 
be done by both adapting the single-particle functions (pj 
as well as by mixing different configurations $j or — in the 
symmetrized version — |n), depending on how close in en- 
ergy they are. As g 00, this culminates in the system's 
fermionization, where there is a one-to-one correspondence 
to fermionic states, n'^ — 0, 1. In particular, the density 
profile p = J2a \4'a\^ becomes broader, with a length scale 
of order \/2N ifJsll . while the strongly correlated nature is 
captured in the fermionic two-body density ^2(2^1, 0:2) = 

/2, which naturally vanishes 



p{xi)p{x2) - \pi{xi,X2) 

at points of collision. 



Harmonic trap (h = 0) 

Here the nature of the fermionization transition has been ex- 
plored extensively for the ground state ll20ll22ll23ll . As for the 
excited states, a look at Fig.|4]unveils that essentially the same 
mechanisms are at work, as exemplified on the one-body den- 
sity p{x) for different states to. The non-interacting density 
profiles have a Gaussian envelope. This may be seen in the 
plot for g — 0.2, the characteristic shape for the degenerate 
states m = 2, 3 stemming from the fact that the interaction 
Hi selects a linear combination so as to be diagonal within 
that subspace. Upon increasing g, the density is being flat- 
tened, reflecting the atoms' repelling one another. Eventually, 
a fermionized state is reached, featuring characteristic humps 
in the density. As in the ground-state case, these signify lo- 
calization; in other words, it is more likely to find one atom 
at discrete spots Xi. However, here the fermionization pattern 
eludes an obvious interpretation, since these are excited states. 
In particular, now the number of humps need not equal N, as 
can be seen for m = 1. 

A look behind the scene is offered by the two-body correla- 
tion function p2 (Fig.|5]l, which includes p — j dx2P2{', X2) 
by averaging over the second atom. It illustrates nicely 
how the interaction imprints a correlation hole at {xi — X2}, 
which relates to the washed-out profile in Fig. ID A complex 
fragmentation of the {xi , X2) plane can be witnessed as we go 
to larger g, which is different from the very obvious checker- 
board pattern of the ground-state case 1 22] . The latter one pro- 
vided a simple interpretation, namely that the atoms are evenly 
distributed at discrete positions over the trap (up to a Gaussian 
density modulation), but with zero probability of finding two 
atoms at the same spot. Here the atoms are apparently more 
localized in the center. On top of that, if one atoms is fixed at 
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Figure 4: Density profiles of A'^ = 4 bosons in a harmonic trap {h 
0) for the excited states m — 1,2,3 (from top to bottom). 



some xi, one cannot unconditionally ascribe definite positions 
for the — 1 remaining particles as before. 



Double well (h — 5) 

For large but finite barrier heights, the lowest excitations at 
(7 = will be formed by the two-mode vectors {tiq, N — n'^). 
All of these will exhibit similar density profiles since p{x) 
only differs significantly near the trap's center; specifically 
p(0) = no \(j)o{0)\^. This can be seen in Fig.|6l which sum- 
marizes the evolution of the lowest excited states' densities for 
= 4. As the perturbation Hj is turned on, different neigh- 
boring states \n) (of equal symmetry) will be admixed, and 
the profiles are adjusted accordingly. Eventually, they saturate 
in the fermionic limit, featuring a typically broadened shape. 
Note that, for intermediate interactions (e.g., g = 1.3), the 



Figure 5: Two-particle density p2{xi,X2) for N = 4 bosons in a 
harmonic trap. Top row: excited state m = 1, . . . , bottom row: 
m — 3; shown are the interaction strengths g = 0.2, 2.2, 15 from 
left to right. 



aforementioned quasi-degenerate states to = 3, 4 (cf. Fig. [3]) 
indeed reveal an almost identical profile. 



The same goes for the two-body density displayed in 
Fig- 13 for these two states (to = 3,4). Absent any interac- 
tions, \E'4 stands out as it has pronounced off-diagonal peaks 
P2{x,—x) = p2{x,x) = 101 (x) I"*, in contrast to ^'3. Al- 
ready for small g, these are washed out due to admixing of 
neighboring states. For intermediate couplings, g — 1.3, 
both densities become virtually indistinguishable. Eventually, 
also here the fermionization limit is reached, and the densi- 
ties can be discriminated again. That said, it should once 
more be stressed that the respective states are by no means 
similar. Rather, they are described roughly through the or- 
thogonal subspaces spanned by {|3, 1), |1,3)} (m = 3) and 
{|2,2), |0,4)} (to = 4). Figure [1] sheds light on this aspect 
by laying out the evolution of the natural populations ?ia(<?)- 
It is similar for to = 3,4 but not identical. The residual 
weights na>2 are slightly separated from the dominant ones 
(a — 0, 1). However, as evidenced for the ground state |.23,], 
they cannot be neglected because they accumulate densely on 
a logarithmic scale, mirroring the extreme one-particle corre- 
lations imprinted in the course of fermionization. 
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Figure 6: Density profiles of A'^ = 4 bosons in a double well (h — 5) 
for the lowest excited states m — 1, . . . , 4 (from top to bottom). 



C. Crossover from single to double well 

We have come a long way studying in depth the spectral 
properties of a single and a double well. As opposed to the 
ground-state case |22], the link between the two is far from 
obvious. In the harmonic trap, the fermionization transition 
was fairly tame, while in the presence of a fixed barrier h = 5, 
there not only seemed to be a strikingly different level struc- 
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Figure 7: Two-particle density p2{xi,X2) for A'^ = 4 bosons in a 
double-well trap {h — 5). Left: excited state m — 3, right: m — 4; 
shown are the interaction strengths g = 0.2, 1.3, 25 from top to 
bottom. The densities for both states are practically indistinguishable 
for mediate g ~ 1. 
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Figure 9: Evolution of the natural orbitals <^a as /i oo for the case 
N = 3 (g = 0.2). Top: The first symmetric orbital (po is notched at 
a; = 0. Bottom: The antisymmetric one (^i) is barely altered. 
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Figure 10: Crossover of the lowest energies Emih) with varying 
barrier strength hfoi N = ?< bosons at interaction strengths g = 0.2 
(a); <; = 15 (b). The line styles are assigned so as to distinguish the 
different level groups at /i = 0. 



ture to begin with, but also the onset of a zoo of crossings and 
quasi-degeneracies. On that score, it would be desirable to get 
an understanding of the crossover from a single to a double 
well. To this end, we will again borrow some inspiration from 
the simple model of a point-split trap h5{x) IstIi . 

Fkst consider the borderline case g = Q. Then the one- 
particle occupations n are conserved for any parameter h, so 
we can assume number states \n) as eigenstates (glossing over 
the fact that, of course, in the instance of degeneracy, the in- 
teraction operator will pick out a unitarily transformed basis 
diagonal in Hi). Let us start with the harmonic trap {h = 0), 
where the spectrum is arranged in steps of Ao = 1 according 
to En = J2a "■a('^ + 5) ^i^^ '^^c paiticles are distributed over 
the oscillator orbitals (pa- Now let us switch on a central bar- 
rier h > peaked at a; = 0. Then each even orbital a £ 2N 
will be notched at x = 0, until its density \4>a\^ will equal that 
of the next, odd orbital 0a+i- Figure |9]gives an illustration of 
this by displaying the natural orbitals 4>q/i. Along that line, 
the energies will evolve continuously from ea to ea+i = a + f . 
On the other hand, granted that the barrier is supported ex- 
clusively at X = 0, the odd orbitals themselves will remain 
completely untouched. Hence, in the limit h oo, we would 
end up with a doubly degenerate single-particle spectrum (or, 
more realistically, a level gap A/j <C 1), which readily trans- 
lates to a shift of AEn = J2ae2N''^'a X 1 =: riovcn' de- 
pending on how many even orbitals were populated to be- 
gin with. Altogether, as the barrier h is run up, the spectrum 
{N/2, N/2 + 1, . . . } at /i = is expected to transform into 



one with a lowest cluster of 1 + (quasi-)degenerate levels 
pertaining to { | tiq , TV — tiq) } at energies E ^ 3N/ 2, followed 
by another one at 3N/2 + 2. 

It goes without saying that a realistic reasoning should take 
into account the finite barrier width (w = 0.5), but the above 
toy model provides us with a rough picture to understand the 
crossover computed for g = 0.2 in Fig. [TOl' a). Note that the 
sketched metamorphosis inevitably brings about crossings be- 
tween different levels as h oo since, for instance, |0, N) 
is barely altered while \nQ,0,N — Uq) is shifted by about 
AE r-. N. 

The above approach may be readily extended to the 
fermionization limit. All we need to do is construct auxil- 
iary fermion states {|n) \ n'^ — 0, 1} and apply the same 
machinery. However, a look at Fig. [TOl b) (g = 15) makes 
clear that the rearrangement of the levels is not as wild as 
as in the non-interacting case. That is simply because the 
'fermions' can only occupy a level once; hence at — > oo 
the lowest group is made up of one or two states only (for 
even/odd numbers, respectively), followed by a cluster of 
four levels regardless of the atom number. You might notice 
that the second band emerging as ft, ^ oo is not perfectly 
bunched at E{h = 10) ~ 11, but really has a runaway at 
E{h — 10) ~ 10.7. This can be traced back to the inclusion 
of a higher orbital 04 in the fermionic state: in such higher re- 
gions, the spectrum ceases to be perfectly doublet-like, foiling 
our previous considerations. 

For intermediate values of g, in turn, one cannot use the 
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same line of argument since the interaction is in the way of a 
simple single-particle description, and n are no longer good 
quantum numbers. Still, the knowledge of the limiting cases 
highlighted above gives a guideline for the crossover. Gen- 
erally speaking, changing h for any g will affect the energy 
via 

i.e, the coarse-grained density p = p * S^^ about the center 
will be reduced so as to minimize the energy costs. This will 
determine the fate of each state when changing over from a 
single to a double well, thus completing our picture of the 
lowest excitations in double-well traps. 

V. CONCLUSIONS AND OUTLOOK 

We have examined the lowest excitations of < 5 bosons 
in harmonic and double-well traps, based on the numerically 
exact Multi-Configuration Time-Dependent Hartree method. 
The key aspect has been the spectral evolution from the 
weakly to the strongly interacting limit, this way extending 
our previous analyses of the ground state [22, 23] to the low- 
est excitations. Moreover, we have illuminated the crossover 
from a single well to a pronounced double well. 

In the case of a purely harmonic trap, the initially equidis- 
tant and degenerate level structure is lifted as interactions are 
introduced, which distinguish between different states of the 
relative motion. In the fermionization limit of ultrastrong re- 
pulsion, a harmonic spectrum is recovered asymptotically. In 
a double well, the non-interacting spectrum has a lowest band 
composed of + 1 states formed from the (anti-)symmetric 
single-particle orbital, well separated from the next upper 
band. Here, the effect of interactions consists in a complex 
rearrangement of the levels, dominated by level repulsion in 
the perturbative regime. Moreover, some lines virtually ad- 
here to one another as interactions are switched on, despite 
their being very different in character. In the fermionization 
limit, we end up with a lowest group made up of the ground 



state (even atom numbers) plus the first excited state (odd), 
followed by a cluster of four levels for any N , that washed 
out due to the non-doublet nature of the higher-lying orbitals. 

In order to get a better understanding, we have also ana- 
lyzed the underlying eigenvectors '^„i- At bottom, the same 
mechanism responsible for the ground-state fermionization 
could be identified. Stronger interactions imprint a two-body 
correlation hole, signifying a reduced probability of find- 
ing two particles at the same position, and eventually lead 
to localization. This becomes visible in the density profiles, 
which evolve from a Gaussian envelope to a significantly flat- 
ter shape. However, the excited states elude an intuitive inter- 
pretation applicable to the ground state. 

Finally, we have cast a light on just how the spectrum is 
reorganized when splitting the purely harmonic trap into two 
fragments. To this end, we considered the deformation of 
the single-particle orbitals as a central barrier is run up. This 
leads to very different energy shifts depending on the overall 
population of even orbitals or, generally, the average density 
about the trap's center 

With these systematic investigations, we have comple- 
mented the extensive work on the ground state. Numerically 
delicate as it is, our study has been limited to the lowest exci- 
tations and also to at most five atoms with an eye toward com- 
puting time. On the other hand, we hope it will also contribute 
to the understanding of dynamical but also thermal proper- 
ties. In this light, an obvious extension would be to study 
time-dependent phenomena. Here double-well systems have 
proven to be a fruitful model for various phenomena. 
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